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ABSTRACT 

We derive the radii of 16 brown dwarfs observed by AKARI using their par- 
allaxes and the ratios of observed to model fluxes. We find that the brown dwarf 
radius ranges between 0.64-1.13 Rj with an average radius of 0.83 Rj. We find 
a trend in the relation between radii and Tefj; the radius is at a minimum at 
Teff~1600 K, which corresponds to the spectral types of mid- to late-L. The re- 
sult is interpreted by a combination of radius-mass and radius-age relations that 
are theoretically expected for brown dwarfs older than 10^ yr. 

Subject headings: brown dwarfs - stars: low-mass - stars: radius 
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Introduction 



Brown dwarfs are generally defined as the objects that are too light to maintain 
hydrogen fusion in their cores. For the solar metallicity, the upper limit of the brow n 



dwarf mass, i.e., the hyd rogen-burning main sequence edge mass, is 0.08 Mq (JKumar 



Hayashi fc Nakano 



1963 



19631 ). The lower edge of mass range of brown dwarfs overlaps with 



mass of planets. They are distinguished by their formation process instead of mass. Brown 
dwarfs are born in the interstellar medium through processes similar to that form stars. 
Objects that are born in protoplanetary disks should be referred to as planets. 

The radii of brown dwarfs and planets in a context of their evolution are interesting 
subject. Radius of an exoplanet can be derived from photometric transit observations. On 
the other hand, it is difficult to determine a radius of a brown dwarf observationally, and 
that has not been discu ssed so far except for the objects in binary system. According to 



the theoretical study by 



Burrows et al. 



(120011 ) some brown dwarfs are actually smaller than 
Jupiter. The radii of brown dwarfs older than 10^ yrs remain roughly constant at 0.7-1.1 Rj 
{Rj: Jupiter's radius) for a broad mass range from 0.3 to 70 Mj. Interestingly, more 
massive dwarfs have smaller radii. This fact is a consequence of the competition between 
two equations of state; Coulomb and electron degeneracy effects. Since the objects held 
by Coulomb effect would only set a fixed density, the relation of radius and mass follows 
r oc M^/^. On the other hand, electron degeneracy results in r oc M~^/^ . Competition 
between two effe cts leads to the roughly constant radius-mass relation at > 10^ yrs. 
Burgasserl (J200ll ) evaluated the expected radius of an observed brown dwarf by perfo rming 



a Monte Carlo simulation adopting the evolutionary models by 



Burrows et al. 



(119971 ) with 



assumption of a constant birth rate and mass function -^ ex M . He found that the most 
likely radius is 0.9 ±0.1 Rj. 



Effective temperature (Tes) of a brown dwarf has often been estimated empirically 
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by assumi n g its radius. IVrba et al.l fl2004f ) used a constant radius of r = 0.9 R,j following 



Burgasserl (l200l[ ) to estimate the e mpirical Tpff of 40 
from the equation L = 47r/2^o"Tjj (JDrilling Sz Landolt 



and T dwarfs. They estimated T^s 



2000l ). where luminosity is derived 



from ii'-band flux by bolometric corrections. They reported that their derived Tcs for the 
early-L dwarfs, about 2400-2500 K, are warmer by about 200-300 K than some earlier 
estimates by fitting the spectral energy distributions with the synthetic spectra (e.g., 



Leggett et al 



200 ll ). They argued that it is caused by different photometry database and 



slightly different assumptions of brown dwarf radii. 



Yamamura et al. (2010) analyzed 2.5-5.0 /im 



Japanese infrared astronomical satellite AKARI ( Murakami et al. 



spectra of six brown dwarfs taken by the 



Tsuji 



2002 



20071 ). They found that 



20051 ) fitti ng to the observec . 



VrbaetaL 



fl2004h . 



Teff derived by the Unified Cloudy Model (UCM; 

spectra are lower by typically 200 K from the empirical Teg estimated by 

They also reported that the radii derived from the AKARI observations have an average 

value of 0.81 Rj, and distr ibute in a wider range between 0.68-1.18 Rj than that expected 



from 



Burrows et al. 



20011 ). They argued that the radii of brown dwarfs should not be 



represented by a single mean value. 

In this paper, we discuss the radii of an extended sample of brown dwarfs determined 
by the AKARI observations in detail. First, we introduce our objects in Section O and 
outline our fitting procedure briefiy in Section |3]q We describe a method of deriving radius 
of brown dwarfs and uncertainty of the radius in Section HI We present the resulted radii 
of 16 observed brown dwarfs in Section 0, and discuss the relation between radius and Tcfj 
and their evolutionary status in Section [6l 



^The details of da t a red uction, calibration and fitting process are described in 



Sorahana &: Yamamura! (120121 ). 
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2. The AKARI Sample 



Twenty seven brown dwarfs including 16 L dwarfs and 11 T dwarfs were observed by 
AKARI. Our targets selected by their expected fluxes to be bright enough for Infrared 
Camera (IRC) onboard the AKARI to provide high-quality spectra within the reasonable 
amount of exposures and their spectral types to sample various types from L to T. We 
obtained good quality spectra of 16 brown dwarfs included 11 L and 5 T dwarfs (Tabled]). 

These 16 objects are nearby and bright, thus they are generally well studied. 
There are three binaries in this sample, GJ lOOlB, 2MASS J1523+3014 and Gl 570D. 
2MASS J1523+3014 is as known as Gl 584C and is a companion of Gl 584AB, which 



is a G dwarf doub 



(JKirkpatrick et al. 



e. 2MASS J1523+3014 is widely (194") separated from Gl 584AB 



2000l ). Gl 570D is a companion of Gl 5 70 ABC triple system, and is 



also located 258" from Gl 570ABC (JBurgasser et al. 



2000l ). In our observation, the target 



source was placed in the 1x1 arcmin^ aperture of the AKARI /IRC instrument. Therefore, 
2MASS J1523+3014 and Gl 570D were observed without confusion from their primary stars. 
On the other hand, G J lOOlB is located only 18.2" from the primary M dwarf GJ lOOlA 



( IGoldman et al 



19991 ) and the spectrum of GJ 100 IB was contaminated by a shoulder of 



an intense signal from GJ lOOl A. We evaluated and subtracte d the signa' 



at the position of GJ lOOlB (see 



Sorahana fc Yamamura 



2012 



for detail). 



from GJ lOOlA 



Golimowski et al. 



( 12004| ) found that GJ lOOlB also has a companion GJ lOOlC separated by 0.087". The two 
dwarfs are of the similar spectral types. To derive the radius of GJ lOOlB (or C; hereafter 
GJ lOOlB), we assume that the luminosity of each dwarf is a half of the observed one. All 
other objects in our sample are believed to be single sources. 
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Table 1. Sixteen Brown Dwarfs observed by AKARI 



Object Name 



Sp. Type T(,ff[K] Parallax(error) [mas] Binary References 



2MASS J14392836+1929149 
2MASS J00361617+1821104 
2MASS J22244381-0158521 
GJ lOOlB 

SDSS J144600.60+002452.0 
SDSS J053951.99-005902.0 
2MASS J15074769-1627386 
2MASS J08251968+2115521 
2MASS J16322911+1904407 
2MASS J15232263+3014562 
SDSS J083008. 12+482847.4 
SDSS J125453.90-012247.4 
SIMP J013656. 5+093347.3 
2MASS J05591914-1404488 
Gliese 570D 
2MASS J04151954-0935066 



LI 


2100 


69.6(0.5 


No 


1, a 


L4 


2000 


114.2(0.8 


No 


2, a 


L4.5 


1800 


85.0(1.5 


No 


1, b 


L5 


1800 


76.9(4.0 


Yes 


1, c 


L5 


1800 


45.5(3.3 


No 


2, b 


L5 


1800 


76.1(2.2 


No 


2, b 


L5 


1800 


136.4(0.6 


No 


1, a 


L6 


1500 


95.6(1.8 


No 


2, b 


L7.5 


1500 


63.6(3.3 


No 


2, b 


L8 


1600 


57.3(3.3 


Yes 


2, b 


L9 


1600 


76.4(3.4 


No 


2, b 


T2 


1400 


75.7(2.9 


No 


2, b 


T2.5 


1400 


6.4(0.3 


No 


3, d 


T4.5 


1200 


95.5(1.4 


No 


2, b 


T8 


700 


169.3(1.7 


Yes 


2, a 


T8 


700 


174.3(2.8 


No 


3, b 



Note. 



Re ferenc e of spectral type (1) iKirkpatrick et al.l \200& \. (2) ICeballe et al.l l|2002h . (3) 



Burgasser et al. 



1 20061 ') 



The pa rallaxes are referred from (a) 
and (d) 



Dahn et al. 



Artigau et al 



1 200611 



1200^), (b) 



Vrba et al. 



1200^), (c) 



Henry et al 



1 20061) . 



The number given for SIMP J0136+0933 is a photometric distance [pc] estimated by comparing the 
spectral energy distribution with known brown dwarfs of similar spectral types. 
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3. Best Fit Model Derived from AKARI Near-Infrared Spectra 

We derived physical parameters of brown dwarfs, namely effective temperature Tes, 
surface gravity log g and critical temperature T^r by model fitting with UCM. Tq. is 
a temperature below which the dust disappears by sedimentation or other unknown 
mechanism, and given as an additional parameter in UCM, i.e., the dust would exist in the 
layer with Tcr< T < Tcond- ^cr is not predictable by any physical theory at present and is 
required to be determined from observations empirically. We mainly use AKARI spectra 
in the range of 2.5-4.15 /im (not to 5.0 /im because the curren t model doe s not explain th e 



observ ed spectra beyond 4.15 /im. See 



Yamamura et al 



20101 ). We follow 



Gushing et al 



( 120081 ) and evaluate the goodness of the model fitting by the statistic G^ defined as 



Gu 



n 

— E 

— m ' ^ 



W,; 



n — 771 



i=l 



fi — CkFf 



kA 



CTi 



(1) 



where n is the number of data points; rra is degree of freedom (this case m = 3); Wj is the 
weight for the i-th wavelength points (we give the equal weight as cjj = 1 for all data points 
because of no bias within each observed spectrum); fi and F^ j are the flux densities of the 
observed data and k-th model, respectively; ctj are the errors in the observed flux densities 
and Ck is the scaling factor given by 



C, 



Y. uJifiFk,i/at 



(2) 



Gk is equivalent to reduced x^? since we adopt cuj = 1 in our analysis. It is difficult to 
determine a unique best-flt model for each object because of relatively large error associated 
with the AKARI spectra. Thus we also use the shorter wavelength spectra (IRTF/SpeX 
an d UKIRT/CGS4 data) t o comp lete our analysis. Details of fltting evaluation are described 



m 



Sorahana fc Yamamural (120121 ). We show an example of model fitting in Figured! This 



object, 2MASS J2224-0158 is a relatively warm object of the spectral type L4.5. 



>. 



> 



1.5 



1 



ao 0.5 

O 







2MASS J 2224-0158 (L4.5) 





reff/log g/rcr = 1800K/5.5/1700K 







1 



2 3 4 

wavelength [|Lim] 



Fig. 1.— The observed and the best fit model spectra for 2MASS J2224-0158 (L4.5). The 
black line is the observed data and the red line is the best-fit model spectra. The spectrum 
between 2.5 and 5.0 /xm is taken by AKARI and that between 1.0 and 2.5 /im is taken by 
IRTF/SpeX. The parameters of the best fit model are given in the legend of the figure. 
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Radius of Brown Dwarfs 



4.1. Derivation of the Radius 

The ratio of model flux F^ [erg/s/cm^/Hz], which is for a unit radius at a unit distance, 
to observed flux f^ [mJy] is written as 



log(F,//, 



-2 Iog(r/L')- 26.497, 



(3) 



where r is the radius and D is the distance. D of each object has been calculated from 
its parallax given in the 4th-column of Table [H except for SIMP J0136+0933 for which 
the value is a spectrophotometric estimate. Distance of SIMP J0136+ 0933 is estimated by 



comparing the flux with known brown dwarfs of similar spectral tvp es flArtigau et al 



he m odel flux Fy is from the UCM best flt model spectra given by 



2006L 



Sorahana &: Yamamura 



( 120121 ). Equation ^ is written with C of the best flt model as, 

C /r_Y 



(4) 



4.2. Uncertainty of the Radius 

The radius in this analysis relies on the distance D and the flux scaling factor C 
deflned in Equation ([2]). The uncertainty of D is estimated from the parallax error shown 
in Table [U The maximum error of D is ~7 %. Some brown dwarfs have multiple parallax 
measurements and the results are consistent to each other within the error. The uncertainty 
of C depends on the absolute flux calibration of AKAR I spectra and the go odness of the 



model fitting. The former is evaluated to be about 10 % (lOhvama et al. 



validated with L' photometry data in 



Sorahana fc Yamamura 



20071 ) , and was also 



(120121 ). The latter error is 



different for each object. Model spectra do not yet perfectly reproduce the observed spectra 
because of for example an incomplete CH4 line list and unknown properties of dust in the 
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photosphere not yet incorporate d into the model (lYainamura et al 



to find a unique best-fit model. 



Sorahana fc Yamamural (120121 ) discussed the uncertainty 



20101 ) ■ Hence it is hard 



of the best-fit model parameters. The uncertainty should be no better than a half of the 
grid spacing (100 K for Tcs and Tcr and 0.5 dex for log g). To estimate the uncertainty 
we change each of T^s, log g, and T^ by one grid from the best-fit value, and search for 
the "restricted best" model by changing other two parameters following the same manner 
through fitting evaluation. If we do not find any models satisfying Gmin "^ Gk < Gmin + 1 
(here, Gmin is taken from all parameter space in model fitting with AKARI data only), 
then the uncertainty of the parameter should be smaller than the grid spacing. In order to 
estimate the uncertainty of C, we evaluate the ratio of G' for the "restricted best" model to 
G for the best model. The maximum deviation of G' /G from 1.0 is ~30 %. The resulted 
overall error of derived radius is listed in Table [21 Uncertainty of radius ranges between 
5 % and 16 %. The uncertainty is the largest for 2MASS J0559-1404 because of the large 
G error. 

5. Results 



The result is listed in Tabled The average radius of 16 brown dwarf samples is 0.83 Rj 
wit h a standard dev iation of 0.14 Rj. The average is slightly smaller than the value given 
by iBurgasserl (1200 ll ). but is consistent within the uncertainty. The derived radii ranges 



from 0.64 to 1.13 Rj. Figure El shows the relation between Tpff and radii at d i fferent masses 



Burrows et al. 



(1200 ih . and 



from 11 Mj to 0.085 Mq predicted by evolutionary model of 

our results are over-plotted with uncertainties. Our results are consistent with theoretical 

predictions for early-L and T dwarfs. 

From this figure we can discuss mass and age of the objects. Both masses and ages of 
2MASS J0036+1821 and 2MASS J2224-0158 are similar to each other, though their Teg 
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are different. Ages of 2MASS J1439+1929 and 2MASS J0559-1404 are close (relatively 
younger than other objects) to each other, but their masses are quite different. Likewise, 
SDSS J1254-0122 and 2MASS J0415-0935 are in similar ages, but with different masses. 
On the other hand, masses of 2MASS J0559-1404 and Gl 570D are close to each other, but 
ages are not; 2MASS J0559-14 4 is y ounger than that of Gl 570D. The age of Gl 570D is 



estimated by 



Kirkpatrick et al. 



(120011 ) as 2-10 Gyr, and is consistent with our result within 



its uncertainty. 

Our sample includes objects with very small radii. Radii of GJ lOOlB, 
2MASS J1523+3014, SDSS J0830+4828, SDSS 1446+0024 and 2MASS J1507-1624 
are 0.64, 0.65, 0.66, 0.74 and 0.79, respectively. They are out of range of the theoretical 
predictions. Interestingly those small radii objects are all around Tcfi~1600 K, i.e., the radii 
of our brown dwarfs hold minimum at mid- to late-L, though we apply some assumptions 
for GJ lOOlB through the data analysis. 

6. Discussion 

We find that the radii of mid- to late-L objects with Tcfr~1600 K are smaller than those 
of theoretical predictions in absolute scale, and those of early-L and T dwarfs in relative 
scale. In this section, we discuss these results. 

6.1. Validity of Our Result 

First, we verify the absolute values of our resultant radii. A fair comparison with the 
theoretical prediction needs a mass and an age of each brown dwarf. However, it is difficult 
to measure the accurate mass and age of each object observationally unless the objects are 
members of a cluster or a binary. If the object is unresolved binary for AKARI, their radius 
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should be even smaller. Mass can be evaluated from the surface gravity. However, the 
uncertainty of mass of each object is very large because of the larg e uncertainty in surface 



gravity derived by the model fitting (jSorahana fc Yamamura 



2OI2I). 



We can justify our result in another way. If we assume that the radius of 2MASS J1523- 
3014 is the mean value 0.83 Rj, the flux levels should be 1.6-1.7 times higher. This is 
unreasonably lager than the error in the flux calibration, and we regard that the radii of 
mid- to late-L dwarfs should be as small as 0.7 Rj. It is noted that the spectra of mid- 
to late-L dwarfs are most affected by the dust in the atmosphere. The incompleteness in 
the atmosphere model could draw a systematic error in the flux level. However, we do not 
see any particular jumps in the flux levels of the models along the changes of the model 
parameters, indicating that 60-70 % changes in the flux level should not be accounted by 
the model. 



Recently, iBurrows et al.l ( 120111 ) discuss metallicity dependence of bro wn dwarf ra dii 



2011 



In fact we argue that our AKARI sample may have metallicity variation ( iTsuji et al. 
Sorahana et al. in prep.). However, the expected change of radius is as much as several per 
cent still much smaller than the discrepancy we found here. 

We conclude that the radii derived from our observed data are regarded to be real and 
admit there exists a serious deviation between observation and evolutional model. 

6.2. Radius Inversion 

Our AKARI sample is selected by their spectral types to sample various types from 
L to T, and may be biased to some extent. They locate very close to our sun (5-25 pc), 
thus we expect that they are as old as the sun (~10^'^ yrs). Ages of thre e binaries in our 



sample, GJ lOOlB, 2MASS J1523+3014 and Gl 570D, are determined by 



Kirkpatrick et al 
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( l200ll ) to be 1 Gyr, 1-2.5 Gyr and 2-10 Gyr, respectively. The ages of our objects can also 
be evaluated from Figure [2|, as we discussed in this section. Those of mid- to late-L dwarfs 
are out of the theoretical evolutionary tracks but they are likely older than 10*^ yr. Our 
results indicate older ages for GJ lOOlB and 2MASS J1523-I-3014. While theory allows a 
wide variet y of mass for T dwar fs, all L dwarfs at near solar age are massive objects (See 



Figure 8 of 



Burrows et al.l 



20011 ). In other words, it is highly possible that L dwarfs of our 
sample are more massive than our T dwarfs. On the other hand, Figure |2] tells that early-L 
dwarfs in our sample (e.g., 2MASS J439+1929, 2MASS J0036+1821) are relatively younger 
and less massive than mid- to late-L objects. 

Brown dwarfs shrink slowly as they evolve. As previously noted, the radii of brown 
dwarfs do no t follow a monotonic function of mass, in particular after the age of 10^'^ yrs. 



According to 



Burrows et al. 



(120011 ). the radius of a less massive object is already small at 
its formation and change of radius during its lifetime is also small. On the other hand, the 
radius of a massive object is large when it is young, but become smaller than the lighter 
objects beyond ~ 10^ yrs, because the electron degeneracy effect overcomes the coulomb 
effect. This inversion of radius takes place continuously in the mass range between 0.3 
and 70 Mj. The relation of radius and mass with age is shown in Figure |3l The data are 
provided from Adam Burrows (in priv. comm.). We see a depression on the curves for the 
dwarfs older than 10^ yrs, which indicates the radius inversion by the degeneracy effect. A 
sharp bump at ~0.015 M© is due to a deuterium burning. The horizontal axes of Figure |2] 
and [3] are not simply equivalent, because effective temperature Teg (spectral type) of a 
brown dwarf depends on both mass and age. According to the discussions above, the trend 
in our AKARI sample can be explained by a combination of both effects. Early types that 
the effective temperature ranges between 1500 and 2100 K tend to depend on their age as 
shown in Figure El We show this trend by overlaying a red arrow in Figure 131 On the other 
hand, late-types with Tefj lower than 1500 K is more mass dependent than age. The trend 
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is described by a blue arrow in Figure |3l Our result of the relation of radii against their Tes 
implies that inversion of radius predicted by theory is actually taking place. 

7. Conclusion 

We derive the radii of 16 brown dwarfs observed by AKARI by comparing model 
flux and observed flux at given distances. The resulted radii ranges in 0.64-1.13 Rj. The 



average radius is 0.83 Rj with a standa rd deviation o 
with the theoretical radii calculated by iBurrows et al. 



0.14 R j. Our results are consistent 



( I2OOII ) within the error for early-L 



and late-T dwarfs, however, a large discrepancy is found for some mid- and late-L dwarfs. 
We find that the radii of mid- to late-L dwarfs are smaller than theoretical prediction. We 
verify that our estimates of radius is reasonable, and conclude that there are deviation 
between observation and evolutionary model. We can not find any reason of this deviation 
and leave it to the future studies. We also find that the radii reach a minimum for the mid- 
to late-L dwarfs. Theory predicts that there is an inversion in the radius-mass relation in 
the brown dwarfs older than 10^ yrs. Our results confirm that this theoretical prediction of 
radius inversion actually present. 

We thank to the anonymous referee for critical reading of our article and constructive 
suggestions. This research is based on observations with AKARI, a JAXA project with 
the participation of ESA. We thank to Prof. Takashi Tsuji for his kind permission to 
access the UCM and helpful suggestions. Prof. Adam Burrows kindly provide us the 
numerical data and warm encouragement. We acknowledge JSPS (PI: S. Sorahana) and 
JSPS/KAKENHI(c) No. 22540260 (PI: I. Yamamura). 
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Table 2. The radius of AKARI objects 



Object Name 


Sp. Type 


Radius (error) [Rj] 


Number in Figure |2] 


2MASS J1439+1929 


LI 


1.01(0.05) 


1 


2MASS J0036+1821 


L4 


0.88(0.05) 


2 


2MASS J2224-0158 


L4.5 


0.87(0.05) 


3 


GJ lOOlB 


L5 


0.64(0.05) 


4 


SDSS J1446+0024 


L5 


0.74(0.06) 


5 


SDSS J0539-0059 


L5 


0.82(0.05) 


6 


2MASS J1507-1627 


L5 


0.79(0.06) 


7 


2MASS J0825+2115 


L6 


0.76(0.09) 


8 


2MASS J1632+1904 


L7.5 


0.74(0.10) 


9 


2MASS J1523+3014 


L8 


0.65(0.09) 


10 


SDSS J0830+4828 


L9 


0.66(0.08) 


11 


SDSS J1254-0122 


T2 


0.84(0.05) 


12 


SIMP J0136+0933 


T2.5 


0.80(0.12) 


13 


2MASS J0559-1404 


T4.5 


1.13(0.18) 


14 


Gl 570D 


T8 


1.04(0.05) 


15 


2MASS J0415-0935 


T8 


0.94(0.05) 


16 
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500 1000 1500 2000 

reff[K] 



2500 



Fig. 2. — The radii of AKARI observed brown dwarfs are plotted on theoretical prediction of 
the relation between Teg and radii at different masses from 11 Mj to 0.085 Mq. Red are less 
massive objects with 11, 12, 13 and 15 Mj, blue are massive objects 0.08 and 0.085 Mq and 
green are intermediate mass objects with 0.02, 0.025, 0.03, 0.035, 0.04, 0.045, 0.05, 0.055, 
0.06, 0.065, 0.07 and 0.075 Mq. The numerical data of evolutionary model are provided by 
Adam Burrows (in private communication). The numbers are object index in Table [2l 
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lO^^yr 
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♦ <> 
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0.02 0.04 0.06 0.08 

Mass [M©] 



0.10 



Fig. 3. — Theoretical prediction of the relation between radius and mass of brown dwarfs 
at different ages. The radii of relatively light and young objects are small and rather stable 
throughout their lifetimes. On the other hand, the radii of massive dwarfs > 0.06 Mq are 
large when they are young, but ultimate radii are smaller than those of less massive objects. 
The radii of early-type objects in our sample tend to depend on their age (shown by red 
arrow), while those of late-type objects tend to depend more on their mass than their age 
(shown by blue arrow). The data are provided by Adam Burrows (in private communication). 
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